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d Abstract 

^ We study the evolution of a many-particle system whose wave function 

obeys the iV-body Schrodinger equation under Bose symmetry. The sys- 
tem Hamiltonian describes pairwise particle interactions in the absence of 
an external potential. We derive apriori dispersive estimates that express 
the overall repulsive nature of the particle interactions. These estimates 
^s^j hold for a wide class of two-body interaction potentials which are indepen- 

(f-) dent of the particle number, N. We discuss applications of these estimates 

£N| to the BBGKY hierarchy for reduced density matrices analyzed by Elgart, 

Erdos, Schlein and Yau. 

qq Keywords: Dispersive estimates; interaction Morawetz type (correlation) esti- 

mates; many-body Hamiltonian; many-body Schrodinger equation; quantum 
hydrodynamics; BBGKY hierarchy. 
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1 Introduction 

a 

Quantum mechanics dictates that the dynamics of ./V interacting particles at 
zero temperature be described by a wave function satisfying the time-dependent 
A^-body Schrodinger equation. When N is large this description becomes im- 
practical. A common remedy is to reduce the system degrees of freedom by 
replacing the many-body Schrodinger equation by effective partial differential 
equations for single-particle wave functions. Thus, the many-body interaction 
is approximately replaced by sums of effective one-particle potentials. This ap- 
proach amounts to a "mean-field" approximation, which is expected to be exact 
as N — > oo in an appropriate sense. This formulation yields nonlinear evolution 
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laws in 3 + 1 dimensions which often lead to successful predictions. In particu- 
lar, the dynamics of particles with integer spin (Bosons) has been described by 
nonlinear Schrodinger-type equations [12j [13l EH [301 E] • 

The relation of mean-field approaches to the many-particle Hamiltonian evo- 
lution is an area of active research. In settings with Boson symmetry, fundamen- 
tal questions concern the rigorous derivation of nonlinear Schrodinger-type equa- 
tions from TV-body evolution laws and the formulation of theories that transcend 
mean field [3]. Recently, Elgart, Erdos, Schlein and Yau (EESY) [7j |9j HO] 
derived rigorously mean-field limits for Bosons on the basis of BBKGY-type 
hierarchies for reduced density matrices [2"r?l |2"T] . 

In this paper we study the TV-particle Hamiltonian evolution of a Boson sys- 
tem for large yet finite TV. We derive certain apriori estimates for hydrodynamic 
quantities, and examine some applications of these estimates to BBGKY hierar- 
chies formulated for dilute atomic gases El [10] . Our goal is to understand 
how many-body properties of the Boson system are connected to the mean-field 
limit. 

The model equation for the particle system reads 

- id t i/j = H n i/j , (1.1) 

where ip is the TV-body wave function and Hjy is the Hamiltonian operator. 
This Hjv has the form 

H N = J2 H a + J2 H " b ' a,b = l,...,N, (1.2) 

a a^b 

where a and b label the particles, H a is the part of the Hamiltonian acting solely 
on particle a, and H a t denotes the (pairwise) interaction between the particles 
a and b. For identical particles, H a is described by 

H a = -A a + Vtrap(Xa) (1.3) 

and the interaction H a b has the form 

H a b ■= \V\ nt (x a ~ X b ) , (1.4) 

where x a denotes the vector position of particle a {x a G Kl 3 ), Vt ra p is a trapping 
potential which confines the particles in space, and Vint is a two-body potential. 



In our formulation, ( 1.1 ) is converted exactly to a set of hydrodynamic equa- 
tions. Apriori estimates for finite TV are derived here directly from these hy- 
drodynamic equations. Note that by we have i/j(t) = U(t)ip(0) where 
U(t) = e ttHN , a unitary operator. For large TV, this representation of ip(t) is 
not particularly useful for apriori estimates and testable predictions. 

We proceed to discuss elements of the model. The trapping potential, Vt rap , 
is introduced in experiments in order to keep the particles together. To avoid 
unecessary complications such as the effect of bound states in our analysis, we 
set Vtrap equal to zero, Vtrap = 0. The case with a nonzero external potential, 
Vtrap ^ 0, is the subject of future work. 
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Next, wc comment on the particle interactions, which are crucial for the 
system evolution. In our model, the pairwise interaction, Vint, is considered as 
a known short-range potential. For dilute gases, which are the primary focus of 
this investigation, the particle collisions are sufficiently rare and weak; hence, 
the distance \x a — xj,\ (a ^ b) remains mostly large in an appropriate sense. 
Consequently, each particle is influenced only by gross features of Vint- The 
specific details of this potential should be immaterial for our purposes. In this 
vein, by considering particles that repel each other, it suffices to assume that 

H ab = V{\x a -x b \) ; V>-c, V'{s)<0 VseR+. (1.5) 

Note that we consider only two-body interactions, neglecting effects of three- 
body and higher-order collisions [29 . 

A key ingredient of the particle model is the scattering length I of the inter- 
action potential, Vint- This I encapsulates the gross features of the two-particle 
scattering [8] . Let the zero-energy scattering solution for the class of potentials 
Vint (5?) that decay sufficiently fast at infinity have the form 1 — wq(x) where 
limi^i^oo Wq — 0; then, I is defined by I := lim| x\^oo(\x\ w o)- Accordingly, Vint 
can be expressed in the scaled form [51 [5TJ [22] 

V int {x) = r 2 V 1 {\x\/l) , (1.6) 

where V\ is a fixed potential with scattering length equal to unity. Thus, the 
Hamiltonian Hn has two parameters, namely, the number of particles, N, and 
the scattering length, In the EESY formulation pH 00 EO HH] , it is essential to 
allow for N — > oo by keeping Nl fixed. A parameter that enters the mean-field 



limit of (jOJ is 

g := AttNI . (1.7) 

Thus, setting g = 0(1) implies that the scattering length I scales with A -1 . 
The recent work by EESY [7J [HI GO HH] deserves special attention, because it 



provides a rigorous justification of the mean-field limit of (1.1 1 for Bosons via 
an averaging procedure. In Sec. [4] we outline their approach and derive certain 
apriori estimates. The uniqueness of the solution to the BBGKY hierarchy in 
the mean-field limit is also proved via a different argument by Klaincrman and 
Machedon [T6"| . 

A crucial observation within the EESY framework is that, in order to be able 
to pass rigorously to the limit as N increases, the two-body potential Vint has to 
be appropriately scaled with N. For example, EESY P [8j [9l [TO] posit that Vi nt 



approaches a delta function as N — > oo, which comes from (1.6 1 with I = N^ 1 . 
These interactions cannot be treated as a perturbation since they have an appre- 
ciable effect for sufficiently large N. We note in passing that in a more elaborate 
model Vi n t is replaced by the Fermi pseudopotential [21 [21 [T3J [T7J [TH1 US], an 
operator that consists of a delta function times a space derivative acting on 
tp. The use of this pseudopotential allows for interesting extensions but causes 
mathematical difficulties and lies beyond our present scope. 
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A few remarks on a mean- field limit of (1.1 1 and its variants are in order 



This limit can be derived heuristically under (1.6) with I = 0(N 1 ) by approx- 
imately replacing the wave function ip by the tensor product 



N 



(1- 



o=l 



where <j) a := (p(t,x a ) and J dx \(f>\ 2 = 1. By manipulation of (1.1 1 and ne- 



glect of terms inconsistent with (1.8 1, a cubic nonlinear Schrodinger equation is 



recovered for <f>(t,x). This equation reads |llj 

-id t ^ = H<t> + g\<t ) \ 2 dp- ii{t)<t> , 



(1.9) 



where H = —A + V tTap and (j,(t) = (g/2) f dx |</>| 4 ; cf . [HI E2 E3 ffl [3TJ • This 
form of /i does not appear in the results of Gross [TH [TJ] , Pitaevskii [25] and 
EESY [71 HJH CEO] but is in agreement with the derivation by Wu E01 EE] ■ This 
fi(t) is not essential since it can be absorbed into a phase factor for <f>. 

It is worthwhile mentioning the complementary view on particle dynamics 
that invokes operator-valued distributions [1, 30, 31]; for the periodic case, see 
e.g. [H][l71[T8]. In this context, one defines the Boson annihilation operator 
and its adjoint, \E'*(x), for a particle at point x. Accordingly, with 
the Fermi pseudopotcntial is written in terms of W and W*. The annihilation 
operator eto(i) for the Bose condensate, the macroscopic quantum state in Bose- 
Einstein condensation, is introduced as the average ao(t) = N~ 1/>2 f dx {(j)*^}, 
where (f>(t, x) is the (one-particle) condensate wave function. The evolution 
equation for <fi can be derived by linearization of Hjy in ^ — a , which is consid- 
ered as small in some sense, and enforcement of (1.1 1 under (1.8 1 pMl I3"U1 I5T]. 



This approach offers physical insight and is amenable to extensions, particularly 
the inclusion of higher-order scattering processes [3] [T7J [THJ [53 [30] . However, 
this formalism is less amenable to rigorous treatment. In particular, the intro- 
duction of the Fermi pseudopotential, which is a type of singular distribution, 
renders the analysis especially difficult. 

Of course, tensor products such as (1.8 1 are not approximate solutions in 
any classical sense, 
meaningful. 



So, an issue is how to make this type of approximation 
The BBGKY hierarchy invoked by EESY [3 [TO] , for which 
the Gross-Pitaevskii equation provides a closure condition, offers a meaningful 
scheme for justifying (1.8 1: if the Gross-Pitaevskii equation is satisfied then all 
the equations in the hierarchy are satisfied as N — > oo. 

To understand the system evolution by connecting macroscopic variables 
such as 4> with microscopic quantities for finite N, we resort to apriori disper- 
sive estimates which express the particle repulsions. These estimates are scale 
independent: we derive them for interaction potentials of the form V(|:e — 
omitting the scaling which can be inserted in the end. Our basic underlying 



assumption is given by (1.51 



Our apriori estimates are useful for the passage to the limit N — + oo, but 
also offer insight into the case with finite TV. The techniques employed in the 
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present work originated in [U [23] in connection to the nonlinear Schrodinger 
equation; see also [SJ I2B] ■ A new element of the present problem is the presence 
of a large number of interaction potentials and the fact that the hydrodynamic 
momenta are not conserved by the flow. Our analysis addresses the question in 
what sense the potentials are repulsive and what are the precise implications of 
this repulsive nature. 

The static version of (1.1) is not addressed here. This case was studied by 
Dyson [6j and by Lee, Huang and Yang [17]. A mathematical proof of the Bose- 
Einstein condensation for the time-independent case was provided recently by 
Lieb et al. [H ESB [HI 122] ■ 

The remainder of the paper is organized as follows. In Sec. [2] we derive an 
apriori space-time estimate for the particle density function after the collapse 
(identification) of two particle positions. Boson symmetry is not required for the 
main result of Sec. [2] In Sec.[3]we restrict attention to Bosons and derive another 
estimate for the square of an appropriately averaged (reduced) particle density. 
In Sec. [4] we outline the EESY approach based on the BGGKY hierarchy for 
reduced density matrices [7[ [5] M EH EH; an d derive apriori estimates in the 
context of this formulation. The (Einstein) summation convention for repeated 
indices is employed throughout this paper unless it is noted otherwise. 



2 First Estimate by Coordinate Collapse 

In this section we derive an apriori space-time estimate for the particle density 
of the iV-body system described by ( 1.1 1-( 1.4 1 when N is large and finite. Boson 
symmetry is not required for our main result but simplifies it considerably, as 
shown below. Our motivation is to investigate in what sense the wave function 
%j) disperses with time. The methodology followed here originates from the work 
of Lin and Strauss [23] . 

Specifically, we show that 



E 



dtdX ab p{t,X ab )<N 2 U\\ m U\\ L 2 



(2.1) 



where the particle density function p is defined by 

1 



p(t,xx, . ..x N ) 



■ x N ) 



(2.2) 



The reduced vector X ab entering (2.1l comes from a basic collapse mechanism 
which identifies two of the variables (a?i, . . . xjsr). Hence, X ab is defined by 

X ab := {X) Sa=Sb e M 3(A, - 1} , X := (x u ,x N ) ; (2.3) 

see Sec. [4] for a discussion on the motivation for this coordinate collapse. For 
Boson particles, tp remains invariant under permutations of a? a 's. Thus, the in- 
tegrals entering the sum of (2.1 1 are the same and the estimate (2.1 1 evidently 
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becomes 

m ab := I dtdX ab {p(t,X ab )} < ||^||iyi||V|U= ■ (2.4) 



We proceed to prove (2.1|. Our program consists of the following steps, (i) 
An evolution equation is derived for the particle momentum density, expressing 
the fact that momentum is not conserved by the flow; see (2.12). (ii) This 
equation is contracted with suitable vector fields to yield (2.171, which contains 
"error" terms depending on the nature of particle interactions, (iii) Estimates 
for these error terms are obtained directly, (iv) The contracted equation is 
integrated in space and time to yield (2.1 1. 

We now describe the procedure in detail. The evolution equation (1.1) is 
recast to the form 

id t i(>+(J2H a )'il)+{J2 H *b)i> = (2.5) 

where H a = — A a and H ab = V(\x a — x b \), which implies that our analysis will 
be independent of any particular scaling of the interaction potential. By (1.5 1 
the basic assumption is 

V'(s) < Vs6l+. (2.6) 

The particle coordinates are denoted by x J a where j = 1, 2, 3 and a = 1, 2, ... ,N. 
We consider the Euclidean space with metric 

gf k := 5 jk S ab , j, k = 1, 2, 3 , a, & = 1, 2, . . . , AT , (2.7) 

and use the notation 

Vj : .[', ■ (2.8) 

OXa 

Next, we define the momentum variables and stress tensor. The components 
of the momentum density for the ath particle are defined by 

^:=^(^V J >-V>ViV*) , 3 = 1,2,3, o = l, 2...N, (2.9) 

where ip* is the complex conjugate of ip. The stress tensor is 

<t ■■= (V^)(V fc V) + (ViV*)(V 6 V) ; 3, k = 1, 2, 3, a, b= 1, ... , N . 

(2.10) 

For completeness, we provide the conservation law for the density p, although 
this law is not invoked directly in the proof. The mass conservation statement 
reads 

ft{p}-V?{pj} = 0. (2.11) 

Use of this law will be made at the end of this section and in Sec. |3] 

The evolution equation for the momenta p" stems from differentiation of ( 2.9 1 
with respect to time and use of (1.1 1. The resulting equation reads 

b,bj£a 



G 



where A := ^2 a A a is the 3iV-dimensional Laplacian. Equation ( |2.12 ) is not a 
conservation law because of V on the left-hand side. By defining the weights 
w a b and the 3iV-dimensional vector M a — (M 1 , M 2 , ... ,M N ) according to 



W ab ■= -2V'(\x a - X b \) 



Ma ■= ^2 W o,bT^- 
b,bj£a 



%b\ 



we rewrite ( 2.12| as 

dt{p a } ~ V£{a b feQ } + V a { Ap} + M a p = 0, a = 1, 2, 



(2.13) 



,N . (2.14) 



For algebraic convenience we have written p a := {p\,P2iPt)i the three-vector 
with respect to the group of coordinates for x a ■ Notice that the condition V < 
for repulsive interactions entails w a b = Wb a _ 0. 

Next, we contract the momentum equation, (2.121 or (2.141, with a suitable 
37V-dimensional vector field, Y( a b)- This field consists of N ordered three- vectors 
and is defined by 



Y ( ab) := (0, ... ,0, 



_ft o 

^6 



,0, ^-^,0,...,0). (2.15) 

\Xb ~ X a \ 



This expression means that the sole nonzero three-vectors forming Y^ a b) are 
the ones in the ath and 6th position in the way indicated above. We adopt the 
convention that Y( a b) = ^(ba) and write Y^, b ^ to denote the jth component of the 
vector located in the cth position of Y( a b) (j = 1, 2, 3 and a, b, c = 1, . . . , N). 
The contraction of the ath-particle momentum, p a , with Y a ( c d) equals 



P a j Y a(cd) — - 



Xd , -* Xd X c 
+ Pd ■ 



Xd\ 



\Xd - x c \ 



(2.16) 



where no sum mation over c and d is implied. The contraction of (2.12) with the 
field Y^tcd) (2-151 produces the equation 

= MpTCW " V H^ M )} + V ■ {^ (cd) Ap} 

Tbvi \^a_ (x7 a V j )Ap+ y Wac *lZ^L Y i p. (2.17) 

= 1 



(v fc ^ (cd) )^ 



V^j Y a(cd)) 



We will show that the dispersive nature of (2.171 provides an apriori bound 
for the collapsed density according to (2.1 1. For this purpose, we integr ate ( 2.171 
over the time-slice [0, T] x M. 3N . The last two terms in the first line of (2.171 are 
explicitly integrated out to zero, viz. 

/ didX V b k {a%Yi {cd) } =0 = | dtdX V]{Yi {cd) Ap} . (2.18) 
We will return to the first term, dt{p'jY^ cd ^ , near the end of this proof. 
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Our main task now is to investigate the nature of terms in the second line 
of (2.171. We sum over all contractions with the vector fields Yt cl j\ for c ^ d. 
There are N(N — l)/2 different vector fields. For simplicity we will relabel 
indices by c = a, d — b and a = c. 

First, we focus on the term — ( V?Y^„,n) Ap, which is integrated to furnish 



3 a(cd)> 



c^d 



J AtAX {-(v«^ M) )Ap} , X = (x u x 2 , 



,x N ) 



The identity 



divY( afc ) = V • Y {ab) = 



Xb\ 



(2.19) 



(2.20) 



and the subsequent integration by parts in (2.191 with c— a and d = b yield 
-4 



E 

a=£b 



AtAX 



\x a - x b \ 



Ap 



; E 

a=£b 



dtdX ab {p(t,X ab }} >0 , (2.21) 



where A a (|x Q — Xb| _1 ) = 47r5(af a — Xb) and X ab is defined in (2.3 1. Notice that in 
our derivation thus far we did not have to assume that the particles are Bosons. 
However, it should be borne in mind that in the case with Bosons the results 
are simplified since all these integrals are equal. 
Next, we focus on the last term in (2.171, viz. 



X a X}j 
\x a - X b \ 



j^b | X b x a 

Wb -x a \ 



p(t,X) 



(2.22) 



which by use of (2.131 can be written explicitly as the triple sum 

X c X a Xb X a 



EE 

a^b c 



\X b ~ X a \ 



W cb -, 



Xb 
Xb\ 



\x a ~ X b \ 



p{t,X) . (2.23) 



Terms with c = b or c = a in the first or second sum of (2.23), respectively, add 

^ a, b. 

(2.24) 



up to a manifestly positive term. We turn our attention to terms with c ^ a, b 
Let us pair the term 

Xc X a Xb X a 



Xb - x a 



which comes from the contraction with Y^ ab ^ , with the term 



X b X c 



X b X c 



(2.25) 



which comes from the contraction with Yi 



(6c)- 



The sum of these two terms 



admits a geometric interpretation as follows. Consider the triangle T(x ai Xb, x c ) 
with vertices x a , Xb and x c and denote the corresponding angles a a := a{a), 
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Figure 1: Geometric interpretation of interaction terms from (2.24) and (2.251 
by the triangle T(x a ,Xb,x c ). 



at := a(b) and a c := a(c) as shown in Fig.[T] Evidently, the sum of terms ( 2.24| ) 



and (2.25) equals 

w ac (cos(a a ) + cos(a c )) > , (2.26) 

which has a positive sign. In a similar manner, one can pair the term 

X c %b x a Xb /r> ~— 

Wcb-TZ ' r= , (2.27 

\x c -x b \ \x a -x b \ 
which comes from the contraction with Y^ ab j , with the term 

Wbc&^r ■ , (2-28) 

\Xb X c \ \X a X c \ 

which comes from the contraction with Yi ac ) ■ The representation by the triangle 
of Fig. [l] implies that the sum of the last two terms is 

w c b[cos(ab) + cos(a c )] >0 . (2.29) 
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As a matter of fact, by inspection of the sum of (2.241 and (2.251 we can 
define the potential 



U(x a ,x b ,x c ) := w c 



X\) x a 



w„ 



x c 



X}) x c 



(2.30) 



%a\ X a \ \%a X c \ \*&b X c \ 

Accordingly, the quantity £ from ( |2.22 ) or ( |2.23[ ) is recast to the expression 

a^b c. 

There are three remarks in order: 



Remark 2.1. There are overall N(N — l) 2 terms in the summation (2.311, 
because each contraction produces 2(N — 1) terms and there are N(N — l)/2 
contractions. The number of terms with c = b or c = a, which have positive 
sign, is N(N — 1); thus, the number of the remaining terms is N(N — 1)(N — 2) 
and we separate them into two groups. 



Remark 2.2. Notice the following term, which comes from (2.241 



V'{\x c -x a \) 



Xb 



\Xb 



(2.32) 



This term can be thought of as a three-particle interaction: the particle c in- 
teracts with a and b. If we consider the plane P a ,b that passes through x a and 
is orthogonal to the vector Xb — x a , then in the half space that contains x b the 



expression (2.321 is positive; in contrast, in the complementary half space this 
expression becomes negative. The idea is that (2.32) can be paired with another 



term where the particle a interacts with b and c so that the sum of the two con- 
tributions is positive. This argument gives general qualitative information on 
the particle mutual repulsion. 

Remark 2.3. By replacing the vector fields X in (2.151 with (x a — Xb — d)/\ x a — 
— d\, where d is some fixed vector, we obtain an estimate like the one in (2.4 1 
where the integration is over x a — x*b = d. 



These remarks conclude the investigation of the last term in evolution law ( 2.17 1 
We have asserted that this term is positive. 



The first term in the second line of (2.171 reads 



V k 1 a(cd) 



ka 
T .l'' 



(2.33) 



Evidently, this term is positive, i.e. E > for the given vector fields. 



It remains to collect the results obtained thus far in order to derive (2.4 1. In 



summary, with regard to the second line of (2.171 we showed that the first and 



third terms are positive and the second term yields the positive integral (2.21 1 
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Recall that (2.18 1 hold. By integrating the first term in (2.171 over the time-slice 
[0,T] x R 3N we obtain the estimate 



dt 



JdX{- d t ( P ^ {cd) )} = JdX [(^ (crf) )| t=0 - { P ^ (cd) ) 



Ti 

(2.34) 



Here, we used definition (2.9) for the momentum variables. This statement 
concludes the proof. 

We close this section with the derivation of another estimate expressing the 
dispersion of the particle density p. This estimate concerns the evolution of the 
variance of p. viz. 

^ J dX { P D} > , (2.35) 



where the weight function D(X) is defined by 
D(x 1; x 2 , ... ,x N ) 



a^b 



(2.36) 



This D is intimately connected with the procedure applied above. Indeed, 
the derivative of D reads 



V a £> = 



X a Xfj 



(2.37) 



With the introduction of Y a by 



V^D(x 1 ,x 2 ...x N )=:Y« , (2.38) 

it becomes evident that we actually performed contractions with the vector fields 
Yj 1 for a = 1,2, ... ,N. Estimate (2.351 for p results from the combination 
of (2.361 with the mass conservation law (2.11 1. 

A possible extension of our analysis is the inclusion of a nonzero trapping 
potential, Vtrap- To speculate the related difficulties, we recall that the repulsive 
nature of particle interactions means that the particles tend to disperse so there 
are no bound states. The addition of a trapping potential in the Hamiltonian 
will produce an effective attractive force that keeps the particles together. This 
influence will compete with the dispersive effect described above, creating the 
possibility for trapped states. 



3 Second Estimate by Commutator Operators 

In this section we restrict attention to Bosons and derive an apriori estimate that 
involves a space-time integral of the square of a reduced particle density function. 
We prove this estimate by resorting to vector commutator operators and the 
evolution of a suitably defined action and associated correlation function. 
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Specifically, we show that 

dtdXa-b \[Pa,b{t, Xab)) 2 } < N 2 \\lp\\ H i\\lp\ 



E 

a=jtb , 



3 

L- 



where the averaged, reduced density p a .b is defined by 

Pa,b(t,x a -b) ■= J dX a ,&d£ a+b {p} . 

R3(JV-1) 

In the above, the variables x a ±b denote the center-of-mass coordinates 

1 _ , 

and the reduced coordinates X a ^ are defined by 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



i.e., X a> b stem from (x\, X2, ■ ■ ■ , f/v) with the pair [x a , x b ) being omitted. 

First, we briefly review and comment on the main assumptions and starting 
equations. The many-particle wave function satisfies ( |2.5| , which is recast to 
the equation 

id t tp- A^ + ^^(|a; a -f fc |)^ = (3.5) 

a^b 



where A is the Laplacian in Mr . The key idea in this section is to invoke 
commutator vector operators, which we apply to conservation laws stemming 
from (13 .51) . In addition to condition (2.6 1 on V, we impose Bose symmetry, 



i.e. require that tp remain invariant under permutation of any coordinate pair 

{x a , Xb) • 

tp(t, ... ,X a , ■ ■ • ,Xb, ■■■) = ''Pit, ... ,x b , ... ,x a , ...) . (3.6) 



This property is inherited by the density function defined by (2.2 1, viz. 
p(t, ... ,x a , ... ,x b , ■■■) = p(t, 



j Xb, 



■) 



(3.7) 



Crucial in our considerations are the evolution laws (2.111 and ( 2.12[ ) for 
mass density and momentum. These equations read 



dtp - V b kP k b = , 
d t p« - V\ {o k J + S 3 k S b a ( -Ap)}+ Mfp = , 



(3.8) 
(3.9) 



where j — 1, 2, 3 and a = 1, 2, . . . , N. The relevant quantities, n ame ly, the 
density p, the momenta p", and the stress tensor are defined by (2.2), (2.9l 
and (2.10). In particular, is rewritten as 

af k = p- 1 [p)p\ + (V»(V£p)] . (3.10) 
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It is worthwhile mentioning the conservation of the total energy, which we will 
not use directly here. This law can be expressed as the conservation of the 
integral 



E(t) 



J dX {~tr(<0+X>(l*a 

a3N a ^b 



(3.11) 



This quantity is a constant (dE/dt = 0) and controls the right-hand side of 
(3.1 1 provided that V is bounded below; see (1.5). 

The proof of (3.1 1 is based on the construction of suitable vector commutator 
operators and study of associated evolution laws. Our program consists of 
the following steps, (i) A vector commutator operator, C, is constructed as a 
suitable average over all particles, (ii) An action, Lit), is defined as the inner 
product (C ■ Pit) | p(t)) where P is the 3-ZV-dimensional momentum vector. 

(iii) The evolution of L(t) is described in terms of distinct inner products, Si. 

(iv) Estimates are derived for each Sg. (v) The evolution equation for L(t) is 
integrated to yield (3.1l. 

We proceed to carry out this program. For each particle pair (a, 6), we 
construct the integral operator B a b by 



{P>abf)(Xa-b) ■ = 



l{x' a+b ) l(X^ b ) 



f(X') dX' 



(3.12) 



\ x a-b x a -b\ 

By using these B ab , we construct the 37V- vector commutator operator 

C a , b := (0 ...0, [x a _ h ; B ab ],0 ...0) , (3.13) 

where [A', B] := AB — BA denotes the commutator of A and B, and [x a -b ; B a b] 
in (3.131 is placed in the ath position and acts on functions / according to 



x a - b ; B ab ~\f)(x a - h ) 



Xg-b x a -b 
\x a -b — x a -b\ 



f(X') dX' 



(3.14) 



For our purposes, it is desirable to symmetrize C a]b and replace it by the operator 

C [a , b ) := (0 . . . 0, [x a - b ; Ba], ... 0, [x b -a\ B ba ], . . . 0) = C a>b + C b , a ■ 

(3.15) 

Subsequently, we average over all particles to obtain the 3-/V-vector operator 



C:=^C Ca;6) . (3.16) 

a<b 

We use this C to construct an appropriate action L(t) whose evolution paves 

h this goal in : 



the way to estimate (3.1). With this goal in mind, let 

=1,...JV 



,2,3 



(3.17) 



13 



be the overall momentum vector and C° a denote the components of C . Accord- 
ingly, we consider the action 



L{t) ■= (C-P(t) | p(t)) . 
The time evolution of L(t) is described by 

L(t) = (C-d t P(t) | p{t)) - {P{t) | Cd t p{t)) . 



(3.18) 



(3.19) 



By using the conservation laws (3.8) and (3.9 1 and transferring the operator 



CV on the right-hand side of the inner product, we obtain the equation 

L[t) = S cm + S cv + S ds + S pT . (3.20) 

The terms Si (£ —cm, cv, ds, pr) on the right-hand side express distinct physical 
effects and are defined by 



p- x V>V|p I (VlCQp 



S rv := (p-'ppl | (V b k Ci)p) - U | (V b k Ci)p? 



(divCX-Ap) | p 
(C ■ (Alp) | p 



(3.21) 
(3.22) 
(3.23) 
(3.24) 



The origin of these terms is described as follows. The term S cm is due to the 
compressibility of the fluid described by the conservation laws; S cv is a convec- 
tive term; S,i s is due to dispersion; and S pi signifies a pressure contribution. 



Our next goal is to show that all terms Sg in (3.21 1— ( 3.24 I are positive; 
in particular, Sds reduces to the integral appearing in (3.1 1. First, we de- 
rive alternative expressions for S cm and S cv in terms of appropriate tensor 
products; see ( 3.35[ ) and (3.361 below. To this end, we introduce some addi- 
tional formulas which will be used below. Recall definition (3.16) and write 
(Ci, C 2 , 



c 



Cm) \ the component C a of C reads 



C'a — ^ [x a -b ; B ab \ 

b,b^a 



and has derivative 

V C C Q 



(v Q - b [f a 

The combination of the identity 



E 

b,b^a 



B, 



S k J B ab 



ab])(S C a -S C b 



b-AK- b Bab) 



(3.25) 



(3.26) 



(3.27) 



with the formula 

([<-„ ; vr^-J /)(*.-*) = 

){Xk,a- 



K-b 



R3« 



I, j 
X a-b 

\Xa-b 



: k,a-b) 



f(X') dX' 



(3.28) 



u a-bl 
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yields the expression 

(VrV a _ b ; B ab ]f){x a . b ) = ( r{{x a ^x' a _ b ) f(X') dX' 



(3.29) 



where the kernel r J k is defined by 



r k { x a~b\ x a -b) : ~ \ x a—b x a-b\ 

X [|^a-6 — ^a'-b| 2 <V ~ { x k,a~b ~ x k,a-b)( X a-b ~ X a-b)] ' 

This kernel is positive definite. By the usual tensor convention we have 
r Q b := (r k 3 {x a - b ;x" a _ b )) , j,k = 1,2,3. 

With regard to the momentum vector P we write P — (p 1 , p 2 , . . 
and denote the tensor product of two momentum components by 

p a ®P b =(pp b k ) , j, k = 1,2, 3. 
For a pair (a, 6) of particles we form the linear combinations 



P 



a±b 



1 



(P a ±p l 



a,b=l, 2, 



JV 



(3.30) 

(3.31) 
■,p N ) 

(3.32) 
(3.33) 



which correspond to the center-of-mass coordinates (3.3). Consequently, we find 
the useful relation 

p a - b ® p a - fe = )^(p a ® p a + p b ® p b - p a ® p b - f ® p a ) . (3.34) 
By combining (3.26l~( |3".34[ ), we express S cm and S cv from ( 3.21[ ) and (3.221 as 
S cm = (p- 1 Vp®V P \{VC) P ) , (3.35) 



(^pop) (vc)p\-/p| (yc)p) 



(3.36) 



With (3.231, (3.24), (3.351 and ( 3.36[ ) in mind, we now show that the terms 
Scm, Scv, Sds and S w entering (3.191 are all positive. In particular, the dis- 
persive term Sd s admits a space integral representation directly related to the 
left-hand side of (3.1 1. 

First, we focus on the convective term, S cv . Since the tensor r(x a - b ; x a~b) 
of (3.31 1 is evaluated at two points, we form the associated two-point momentum 



J(X;X') := v^Vp P - ^p/p^ P' , 



(3.37) 



where 



P := P(x a+b , x a - b , X a>h ) , P' = P{x' a+b , x' a _ b , Xi b ) , (3.38) 
p := p(x a+b ,x a - b ,X atb ) , p' = p{x' a+b , x' a _ b , X^ b ) . (3.39) 
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Accordingly, the convective term of ( |3.36| is expressed as the sum 

S cv = 2Y,W ab , 

a^b 



where the integrals W ab are given by 

W ab := J dx a+b dx' a+b dX a , b dX^ b j(j a - b ® J a - fc |r a; 
with the obvious convention 



1-6 



> 



J 



a±b 



1 



( J a ± 



V2 



(3.40) 



(3.41) 



(3.42) 



In (3.41 1, (• | • ) _, denotes the expectation value, or average, over the variables 



(x a -b, x' a _ b ). By (3.401 and (3.411 we conclude that 

S cv > . 



(3.43) 



Second, we concentrate on the compressible term, S cm , from ( 3.35[ ). By 
introducing the gradient operator 



ia±b 



V2 



(V a ± v fc 



S cm is recast to the form 



Scm = 2 (p^ 1 V< " "p <8> V a - V | (r ab )p) 



a^b 

Without further ado, we conclude that 

Scm > 



(3.44) 



(3.45) 



(3.46) 



Third, we show that the pressure term, S pr , is also positive. This case is 
more demanding but crucial since it is connected intimately with the nature of 
particle interactions. Equation (3.241 is recast to the expression 



Spr (M a p) | P )= ]T Q a 



(3.47) 



a,b,c 



where the triple interaction terms Q a ,b,c have the integral representation 

x c x. 



Qa.l 



Sb £a S ^-p(X)p(X') dXdX' . (3.48) 



\Xc X a \ \x b X a {x b % a) 



Recall that w ca is defined in (2.131. To derive the desired estimate, we pair 
appropriately terms participating in the sum of (3.471. So, we pair the terms 

y P (X)p(X>), (3.49) 



X c X a Xfo X a Xfo "T" x a 



\%c %a\ -Eci x b X a\ 
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Xjy X C Xfr ~\~ % c 



(3.50) 



\%c -Ea\ \%b ^b~^~^c\ 

In the spirit of Sec.[2j we interpret the above terms geometrically by considering 



the triangles 

T(x a , x c , x b - x b ' + x'g) =: T(a, c,b-b' + a') , 
T{x a , x*c, x b — x' h + x' c ) =: T(a, c, 6 — b' + c) , 

as shown in Fig. [2} By symmetry relation (3.7 1 for Bosons, we have 

p{. . . x' a . . . x' c . . .) = p{. . . x' c . . . x' a . . .) . 



(3.51) 
(3.52) 



Hence, we can add to (3.491 and (3.501 the contributions from the terms 



X b X a X^ -\- X c 



Cctj \x b X a Xfo -\- X c \ 



7, P( X ) P( X ") . 



X\j X c Xfo ~\- X a 



V a \ \Xfj X c %b~^~%a\ 



71 P(X) P(X") , 



(3.53) 



(3.54) 



where X" denotes the variable resulting after we switch x' a and x' c . Accordingly, 
we consider again the triangles 



T(a, c, b — b' + c) , T(a, c, b — b' + a') . 



(3.55) 



Thus, by the notation of Sec. [2j the sum of the four terms described in ( 3.49| , 
( |3~50| , ( [3~53"| and ( pHilf equals 

w ca {cos a{a, c,b — b' + a') + cos a(c, a, b — b' + c')} 

+w ac {cos a(a, c, b — b' + c') + cosa(c, a, b — b' + c')} ; (3.56) 

see Fig. [2] Evidently, this term is positive. Hence, we conclude that 

S PT > . (3.57) 

We now focus on the dispersive term, 5d s , and invoke the Fourier transform 
with respect to the variables {x a }. Let £ a be the dual variable corresponding 
to the three- vector x a . By virtue of the dual 3 N- vector 

the Fourier transform of the 3iV-dimensional Laplacian reads 



T{-A}=\C +b \ 2 +\ia- 



1 3 I 2 

L — 'a.b 



(3.58) 



Note that formula (3.23) for Sd s involves divC. By (3.12 1-(3.16 1, the divergence 
of C is 

divC = 2j2 B ab ■ (3.59) 
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Figure 2: Geometric interpretation of interaction terms from (3.471 and (3.481 



In view of (3.12), the Fourier transform of the requisite B a b is 

<5(£a+b) <5(2 a ,&) 



\Ha-b\ 



Thus, we readily compute the transform 

^{div(7(-A)} =J2Hia-b) S(C+b) S(E a<l 



(3.60) 



(3.61) 



a^b 



By expressing Sds in terms of Fourier transforms, we conclude that the dispersive 
term is in fact the integral 



S ds =^2 j d %a-b {(Pab) 2 } 
a+b „ 3 



(3.62) 



where the reduced density p is introduced in (3.2 1. 

So far, we proved that the terms entering the right-hand side of (3.191, the 
evolution equation for L(t), are positive. In particular, Sds is described by 
integral (3.62). The final stage of the proof involves integrating (3.201 over a 
time interval [0,T]. Thus, we obtain the equation 

/ dt {S cm + S cv + S ds + S pr } = L(T) - L(0) . (3.63) 
Jo 
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Furthermore, the action L(t) defined in (3.181 is bounded by H^llff 1 HV'llz^- Con- 



sequently, we reach the desired space-time estimate (3.1l. Note that the right- 



hand side of (3.1 1 is bounded by the energy E(t) of the system; cf. (3.11 1. This 
statement concludes the proof. 

We close this section with a few comments. By defining the function 



j-6 



D(X;X') :=J2\xa-b 

a<b 

and the correlation function 

C(t) := JdXdX' {D(X;X') p{t,X)p(t,X')} 



we deduce the relations 



C(t) = L(t) , L{t) > 



(3.64) 



(3.65) 



(3.66) 



Because L(t) is increasing and bounded by the total energy E(t), we conclude 
that 

lim L(t) = L ± , (3.67) 

t — >±oo 



i.e., C(t) in (3.651 is a convex function with a unique minimum, and is asymp- 



totically linear at t = ±oo. 



4 Correlations and BBGKY Hierarchy and Es- 
timate 

This section has two parts. The first part serves a brief review of the description 
of the many-body Hamiltonian evolution for Bosons in terms of a BBGKY 
hierarchy for particle reduced density matrices [TJ [51 El EH 127] ■ I n the 
second part, we derive an estimate for the one-particle marginal 71; see (4.171 



and (4.20) below 



4.1 Review of BBGKY hierarchy for A;-particle marginals 



We start again with the A^-body Schrodinger equation ( |3.5| . The ultimate goal is 
to investigate aspects of this evolution for large N. Following EESY [Tl 151 19"1 ITU] 
we apply the idea that the two-body interaction V is scaled in some way by 
incorporating the notion that pair collisions are relatively rare and weak. One 
particular scaling for V is [71 [HJ GO [TD] 

V(\x a -x b \) = N 2 V!(N\x a ~ x b \) = l{iV 3 Vl(iV|f a -f 6 |)} , (4.1) 



which is equivalent to (1.6 1 with a scattering length I = 0(N 1 ). The scaling 



with N is chosen so that the quantity inside the curly brackets approaches 
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a constant times the delta function as N — * oo. We remind the reader that 
because of Boson symmetry the iV-body wave function ip is invariant under 
permutations of the space variables, i.e. 

1p(t,X!, X 2 , ■■■ ,X N ) = i>(t,& n (i-), f T ( 2 ), • ■ ■ ,Xir(N)) , (4-2) 

where 7r(a) is a permutation of the indices {1,2, ... ,N}. Although we are 
interested in the time evolution of the system, for the sake of convenience we 
will omit the time (t) dependence whenever this is not relevant to the discussion. 

A description of the particle system via a BBGKY hierarchy was proposed by 
Spohn [26JE2] and forms the starting point in the analysis by EESY [71 18} f0| HO] . 
A key element of this methodology is the density matrix jn '■= \' I P) (ip\, which is 
a trace class operator represented by 

j N (t,Xi, ...x N \x[.. .x' N ) := ip(t,Si, . ..x N ) ip*(t,x[ . ..x' N ) . (4.3) 

By virtue of (|3.5|), the evolution equation for jn is 



idtlN + (-A N + A N , ) 7JV + J2 " Va ' b ' hN = , (4.4) 

where An = J2 a ^« denotes the 3iV-dimensional Laplacian, V a b ■— V(\x a — Xb\), 
and the primed indices indicate the coordinates [x[ . . . ^Jv)- 

Next, we describe features of the fc-particle marginals "fk stemming from 
7at. For this purpose, it is convenient to define := [x\, X2, ■ ■ ■ ,a?jv) and 
X' N := {x[, . . . , x' N ). In order to average over some of the coordinates, we also 
define E M US] 

X k = (xi, ... ,x k ) , X N -k = (a?fc + i, ... ,xjy) ; (4.5) 

thus, Xn = (Xk, XN-k) for k — 1, 2, ... ,N. The fc-particle marginals as 
formed via the partial averaging 

7fe (A fc I XI) = f dY N - k {j N {X k ,Y N _ k I X k \Y N - k )} . (4.6) 



Note that because of (4.2 1 it does not matter which variables we average out. 



We now describe the evolution law for 7^ using (4.4 1 and definition (4.6 1 



By (3.3), let V ±& be the gradient operator corresponding to x* a ±b- I n view of 
A a — A a > = V a+a ' • V a _ a ', setting a — a' and integrating over a + a 1 for the 
averaged variables j k yields zero. Define the sets J k ■— {1,2, ... ,fc}, J' k := 
{V, 2', . . . , k'} and J k = {fc + 1, . . . N}, the complement of J k . The above 
averaging also produces zero for the potentials V a b and V a >b< if a, 6 € J k or 
a', 6' G Jf.' when we identify the variables a — a' and b = b' , i.e. we identify the 
coordinates (x k +i, ■■■ ,%n) with (x' k+l , ... ,x' N ). Consequently, the evolution 
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equation for 7^ reads 



ida k - (A k - A' k )jk 



E v »»- E ^ 



x,6£J fc 



',b'e.j' k 



Ik 



a£Jk,a'£j' k b£J k 

xjn(xi .-.x a ... x k ,y k +i ■■■Vb--- Vn\x[ 
= . 



■ x' a ...x' k ,y k+1 ...y b ...y N ) 
(4.7) 



We proceed to write (4.7 1 as a BBGKY-type hierarchy for 7^. Because of 



Boson symmetry, integrations over b £ J k are reduced to 

- fc ) E / dik- fc -idfc + i N** - fc+il) ~ V(K - y k +i\)] 

x~f N (xi ...x a -. .x k ,y k+ i,Y N -k-i\x[ ...x' a . . . x' k ,y k+ i,Y N - k -i) .(i.8) 
With the definitions Hk ■= —Ak + J2 a beJk V a b where a ^ b and 

Cvi+i bk+i] = y dj/fc+i [V(\x a -y k+1 \)-V(\x' a -y k+1 \)]) 

x*y k+ i(X k ,y k+ i\X' k ,y k+ i) , (4.9) 



(4.7 1 takes the form 

idtlk + (H k - H k ,) lk + (N-k)J2 E C vi+i K+i] = ' ( 4 ' 1Q ) 

This set of coupled of equations constitutes a finite BBGKY-type hierarchy for 
7fe where k = 1, 2, . . . , N. Evidently, each 7^ depends on N. 

4.2 Estimate for one-particle density matrix 71 

Thus far, we have essentially reviewed the formulation by EESY [H [10] . In 
the following, we exploit this framework to derive an estimate that involves 71. 
We assume that the interaction potential V\ of (1.6 1 is integrable, V± € L 1 ( 



For finite N, the first equation of hierarchy (4.101 reads 



id tl i - (Ai - A r ) 7l + (N- 1)4£ [ 72 ] = , 



(4.11) 



where 72 = 72 (afi, a?2 | x{, x^) . To motivate the analysis for finite N given below, 
we first consider the simpler case where N — > 00, in which the interaction 
potential approaches a constant times a delta function. As N — > 00, by the 



scaling of (4.1) we have the limit 

V N (-) :=(N-l)N 2 V 1 (N\x a -y k+1 \)^g6(\x a -y k+1 \) , 



(4.12) 
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where g is a fixed positive parameter (<? > 0). In effect, we arrive at the basic 
collapse mechanism that reduces 72 to a function of 6 variables, viz. 



#2 [72] : = 9 {72(^1,^1 I -l2(xi,x[ I x[,x[)} 



Thus, (4.11) becomes 



id t rfi - (Ai - Ai/)7i +^[72] = . 



(4.13) 



(4.14) 



The partial Fourier transform of ( |4.14[ ) has the relatively simple structure 
of a transport law; see (4.18). Next, we derive an estimate for the forcing term 
02 [72] °f (4.141. Notice that the operator B2 acting on 72 is written as 

82(72] = 9 J dY N - 2 ^ip(xi,xi;Y N - 2 ) tp* (xx,x{;Y N - 2 ) - ip(xi,x[;Y N - 2 ) 



By writing 



dY N _ 2 



1 -2 

-\lp(Sl,X2,Y N -2)\ 



(4.15) 
(4.16) 



and denoting 02 [72] by B 2 {xi,x{), we have the estimate 



Ba(t) 



Axi- V 
< 2 dxi-v 



dxi-i 



B 2 (xi + v , xi-i>) 



< 



dfi+i'|Vp(^i, 
dxi-i' 



xi 



dxi p(xi,xi) 
J dxi +v p(xi,x{fj 
dx 1 p(x lt xi) , 



(4.17) 



by the definition of p and p. It should be borne in mind that the density p 
depends on t. 

In order to study the implication of (4.171, let us write Ai — Ay = Vi+i' ■ 
Vi_i' and B 2 (xx,Xi') = 62(^1+1') £i-i')- We proceed to derive an estimate 
involving 7^ By taking the Fourier transform of (4.141 in the x%+y variable 
and denoting the dual variable by v, we arrive at the transport equation 

drfi-2irir-V 1 -i,y 1 =Ba(t,it,2i-v) , (4.18) 

where 71 (i, v, xi-y) denotes the Fourier-transformed function. For the sake of 
some simplicity, we assume that the initial data for 71 are zero, i.e. 71 (t = 0) = 
0. Thus, the solution to (I4.18I) is 



ji(t,v,xi-r) 



02 (s, v, xi-i' — 2irsv) ds 



(4.19) 
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On the basis of (4.19), we derive the estimate 



sup||7i||£3 

t,v 



sup 

t.v 



J ds Jdy {e a ^B 2 {s,y,x^ v -2tt«3)} 



< 



ds 



B2OO 



(4.20) 



The combination of (4.201 with (4.171 and (2.1l provides the desired estimate 
for 71 when N — > 00. 

This anal ysis a p plies directly to the case of finite iV with minor modifica- 



tions. First, (4.18 )— (4.20 1 remain intact under the replacement (N—1)C 



V.2 



B 2 



in (4.11). Second, in view of Remark 2.3 of Sec. [2] and for Vi G L^M 3 ), esti- 
mate (4.17) becomes 



< 



ddVj\r(<i) / da?i /a(a?i , a?x + d) 
ddVN(d) / dxi p(a?i, xj. + d) 



dxi_ : 



< ||^v|| L i sup / dx 1 p(x 1 ,x 1 + d) . (4.21) 

It is natural to ask whether this analysis can be extended to fc-particle 
marginals for k > 1. However, evolution equations (4.101 satisfied by 7^ cannot 
be converted to simple transport laws if k > 1 because of the presence of interac- 
tion potentials. Hence, the procedure applied hitherto is not directly applicable 
to the case with k > 1. The derivation of estimates for "fk>i IS left f° r future 
work. 
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